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UrrtODUCTION 

The  well -publicised  adulteration  of  nonprescription  capsules  with 
poisons  has  called  attention  to  the  need  for  rapid,  noninvasive  and 
nondestructive  aethods  of  screening  over-the-counter  drugs.  In  1082 
potassium  cyanide  appeared  in  capsules  of  Extra  Strength  Tylenol  in  the 
Chicago  area,  and  resulted  in  seven  deaths  (1,2).  Subsequently,  a  nuaber  of 
"copycat"  incidents  occurred,  involving  such  things  as  strychnine,  mercuric 
chloride  and  sodiua  hydroxide  in  capsules,  hydrochloric  acid  in  eyedrops 
and  sodiua  fluoride  in  an  artificial  sweetner  (£~4) •  Concern  about  drug 
taapering  has  continue^  to  grow;  earlier  this  year  cyanide-laced  Tylenol 
capsules  caused  the  death  of  a  woman  in  the  New  York  area  (5,6),  and  even 
wore  recently  Excedrin  capsules  containing  cyanide  resulted  in  two  further 
fatalities  in  the  Seattle  area  (7). 

The  nuaber  of  cases  of  product  taapering  appears  to  be  increasing.  In 
addition  to  analgesics,  other  types  of  capsules,  including  cold  and  allergy 
reaedies  and  an  appetite  suppressant,  have  been  affected  (8,9).  The  cost  to 
Johnson  and  Johnson  alone,  the  maker  of  Tylenol,  is  estiaated  to  be  S150 
million  (10)  for  the  most  recent  recall  of  capsules.  Litigation  concerning 
the  1082  Tylenol  poisoning  is  still  proceeding;  however,  there  is  an 
indication  that  judges  and  juries  alike  are  beginning  to  hold  manufacturers 
responsible  for  taapering,  even  when  there  is  no  evidence  linking  the 
company  to  the  actual  incident  (under  the  doctrine  of  res  tpst  Loquitur, 
i.e.,  the  thing  speaks  for  itself)  (11).  In  such  cases  the  burden  falls 
upon  the  company  to  prove  that  taapering  did  not  occur  in  the  company's 
plant  (11).  It  is  the  objective  of  this  report  to  demonstrate  that  near- 
infrared  reflectance  analysis  (NIRA) ,  combined  with  a  suitable  cluster- 


analysis  algor  it  ha,  can  be  used  to  indicate  that  capsules  (or  similar 
products)  have  or  have  not  been  adulterated. 

Subsequent  to  the  1082  incident  the  Food  and  Drug  Administration  (FDA) 
collected  and  tested  two  million  capsules  of  Tylenol  in  a  search  for 
contaminated  bottles  of  the  drug  (12).  Capsules  were  batched  into  groups  of 
5-10  for  analysis.  While  the  FDA  avoided  much  publicising  of  its  methods, 
several  techniques  have  been  reported  for  various  capsule  analyses, 
including  such  simple  methods  as  inspection  by  visual  appearance  and  odor. 
In  the  Tylenol  case  the  capsules  had  been  grossly  contaminated  with  500  to 
800  of  KCN,  and  the  KCN  consisted  of  fairly  large  crystals  while  the 
analgesic  was  a  powder  of  small  particle  sise  (4,12).  Potassium  cyanide  is 
also  deliquescent,  which  resulted  in  a  readily  identifiable  distortion  and 
discoloration  of  some  of  the  adulterated  capsules.  KCN  can  also  emit  an 
odor  of  bitter  almonds  (13,14).  Nevertheless,  the  possibility  exists  that 
tainted  capsules  will  not  be  spotted  and,  in  any  case,  the  procedure  is  not 
suitable  for  detecting  lower- level  contamination.  In  another  case,  UV 
spectrometry  has  been  used  to  identify  substitution  of  phentersine, 
phenylpropanolamine,  and  caffeine  (15).  Thin-layer  chromatography  (15)  and 
microcrystal  tests  have  also  been  used  to  detect  counterfeit  lonamin 
capsules  (15).  X-ray  spectrometry,  using  grain  inspection  or  clinical 
aammographic  instrumentation,  has  been  employed  to  detect  cyanide  in 
Tylenol  capsules  (12) ;  however,  the  majority  of  these  determinations  were 
performed  using  differential  pulse  polarography  of  cyanide  reduction  at 
about  -0.3  V  vs.  SCE  (12).  Inductively-coupled  plassa  atonic  emission 
spectrometry  (ICP-AES)  of  trace  elements  in  KCN  has  also  been  used  to 
obtain  an  elemental  "fingerprint*  of  the  tainted  capsules  in  an  effort  to 
trace  the  source  of  the  KCN  (4) ■ 


Vitb  the  possible  exception  of  X-ray  Methods,  all  of  the  above  techniques 
require  that  the  capsules  be  opened  and  their  contents  eaptied  for 
analysis.  Clearly,  a  nondestructive  aethod  of  probing  the  contents  of 
suspect  capsules  (or  siailar  products)  directly  through  the  walls  of  the 
container  would  be  desirable  for  rapid  screening  in  large  nuabers.  NIRA,  in 
conjunction  with  appropriate  data  analysis,  is  suitable  for  such  a 
screening  of  over-the-counter  drugs  and  is  siapler  and  less  expensive  to 
iapleaent  than  x-ray  Methods. 

NIRA  is  a  rapid  analytical  technique  that  typically  uses  the  diffuse 
reflectance  of  a  saaple  at  several  wavelengths  to  deteraine  the  saaple's 
coaposition  (16).  Through  a  computerised  Modeling  process  NIRA  is  able  to 
correct  automatically  for  background  and  saaple-aatrix  interferences, 
aaking  ordinarily  difficult  deterainations  seea  routine.  This  Modeling 
process  eaploys  a  'training  set"  of  saaples  to,  in  effect,  "teach*  the 
coaputer  to  recognise  relationships  between  ainute  spectral  features  and 
saaple  coaposition.  Of  course,  the  contents  of  the  training-set  saaples 
aust  have  been  previously  deterained  by  soae  other  aethod. 

The  aodel  developed  in  NIRA  is  coaposed  of  linear  equations  of  the  fora 

n 

Concentration (A)  =  C0  *  l  CjR*  (1) 

i=l 

where  A  is  the  saaple  coaponent  of  interest  (one  equation  is  required  for 
each  coaponent),  n  is  the  nuaber  of  wavelengths,  R*  is  the  reflectance  at 
the  i-th  wavelength,  and  the  C’s  represent  the  constant  factors  deterained 
through  a  Multiple  regression  process.  In  other  words,  the  aodel  gives  the 
saaple  coaposition  froa  a  nuaber  of  linear  equations,  each  of  which 
expresses  a  particular  coaponent  concentration  as  a  weighted  sub  of  the 
reflectances  observed  at  a  nuaber  of  wavelengths. 
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The  instrument  used  in  NIRA  can  be  as  simple  as  a  filter  photoseter  or 
as  complex  as  an  FTIR  (the  former  is  more  common  than  the  latter) .  The 
broad  spectral  features  and  highly  correlated  wavelength  vectors  (R^'s) 
make  only  a  few  filters  (less  than  10)  usually  necessary,  so  NIRA 
instruments  are  relatively  inexpensive.  Little  or  no  sample  preparation  is 
required  in  NIRA,  and  powders  can  be  directly  analysed.  Finally,  near- 
infrared  radiation  penetrates  most  compounds  rather  well  because  the 
absorptions  in  this  spectral  region  are  usually  weak.  These  basic 
characteristics  suggest  that  NIRA  can  be  easily  and  profitably  employed  in 
the  detection  of  tampering. 

However,  there  is  an  obstacle  to  such  a  use  of  NIRA:  it  is  not  possible 
to  predict  what  might  be  placed  in  a  particular  product.  The  modeling 
process  described  above  relies  on  the  availability  of  a  training  set 
composed  of  known  products  and  known  contaminants.  Bven  if  one  could 
assemble  sets  of  all  of  the  products  and  adulterants  that  have  been 
involved  to  date,  there  is  no  guarantee  that  a  new  adulterant  would  not 
appear  in  the  product  tomorrow.  Unfortunately,  using  a  multiple  regression 
model,  any  amount  of  reflectance  at  the  selected  analytical  wavelengths 
generates  some  sort  of  composition  value  regardless  of  the  material 
responsible  for  the  reflection.  In  other  words,  when  a  sample  contains  a 
component  that  is  not  present  in  the  training  set,  erroneous  composition 
values  can  result  without  any  indication  of  the  error. 

One  cure  for  this  problem  would  be  to  find  a  method  of  spotting 
anomalous  samples  based  on  their  near- inf rared  spectra.  The  use  of  such  a 
method  would  allow  different  linear  models  (calibration  equations)  to  be 
applied  in  the  analysis  of  the  components  of  different  samples.  The  problem 
of  assigning  a  particular  spectrum  to  a  particular  linear  model  has  been 


called  the  *false-saaple  problea" ,  and  a  aethod  has  been  proposed  to  solve 
it  (17).  This  aethod,  the  Quantile  BEAST  (Bootstrap  Error-Adjusted  Single- 
saaple  Technique) ,  goes  beyond  a  siaple  qualitative  analysis  of  aixtures  to 
deteraine  whether  a  quantitative  prediction  equation  applies  to  a 

particular  saaple.  This  aethod  can  be  used  to:  (1)  detect  any  taapered 
product  by  deteraining  that  it  is  not  siailar  to  the  previously-analysed 
unadulterated  product,  (2)  qualitatively  identify  the  contaainant  froa  a 
library  of  known  adulterants  in  that  product,  and  (3)  provide  a 

quantitative  indication  of  the  aaount  of  contaainant  present. 

TK01Y 

The  Quantile  BEAST  considers  each  aonitored  wavelength  to  be  a 

diaension  in  hyperspace.  For  exaaple,  a  spectrua  recorded  at  n  wavelengths 
can  be  represented  as  a  single  point  in  n-diaensional  hyperspace, 
translated  froa  the  origin  in  each  diaension  by  an  aaount  that  corresponds 
to  the  aagnitude  of  the  reflectance  observed  at  each  wavelength.  In  this 
scheae,  siailar  spectra  appear  in  siailar  regions  of  hyperspace.  The 
distribution  of  reflectances  on  each  wavelength  axis  provides  a  projection 
of  the  clusters  of  siailar  points  (spectra)  (see  Figure  1).  Valid 

(unadulterated)  saaples  are  defined  as  those  that  fall  inside  the  cluster 
of  training-set  points  when  the  BEAST  is  trained  with  unadulterated  product 
saaples.  False  (taapered)  saaples  are  those  that  fall  outside  of  the  saae 
cluster.  Confidence  liaits  are  set  along  any  linear  coabination  of 
wavelengths  (diaensions)  to  define  the  surface  of  the  cluster  at  a 
specified  level. 


These  confidence  liaits  are  obtained  by  using  a  bootstrap  procedure 
(18)  to  arrive  at  an  estiaate  of  the  real-saaple  distribution  based  on  the 
training-set  distribution.  The  bootstrap  procedure  has  three  basic  steps: 

1.  A  training  set  is  carefully  constructed  froa  real  (in  this  case 
unadulterated)  saaples  in  a  way  that  captures  all  of  the  possible 
saaple  variance.  Fortunately,  NIRA  is  predicated  on  the  construction 
of  such  training  sets,  and  a  good  deal  of  knowledge  has  been 
accuaulated  in  this  area  (19,20).  Spectra  are  then  recorded  to 
produce  training-set  points  in  hyperspace. 

2.  A  randomly  selected  set  of  samples  (the  same  size  as  the  training 
set)  is  chosen  from  the  training  set,  with  replacement  from  the 
training  set,  to  form  a  bootstrap  sample  set.  The  bootstrap  sample 
set  is  analyzed  as  though  it  were  an  actual  sample  set  -  in  the 
tampering  work,  this  analysis  is  based  on  locating  the  center  of  the 
sample  distribution  in  hyperspace. 

3.  The  complete  hyperspace  distribution  of  real-sample  centers  is 
approximated  by  the  bootstrap  distribution  of  centers.  This 
approximation  is  calculated  using  a  Monte  Carlo  integration  of  the 
bootstrap  distribution  while  the  training  set  is  held  fixed  at  the 
values  obtained  in  step  1 . 

After  it  develops  an  estimate  of  the  real-sample  distribution  from  the 
training  set,  the  BEAST  computes  the  center  of  the  real-sample 
distribution.  Vhen  a  new  sample  (suspect  product)  is  analyzed,  its  spectrum 
is  projected  as  a  point  into  this  hyperspace.  A  vector  is  then  formed  in 
hyperspace  between  this  new  spectral  point  and  the  computed  center  point  of 
the  real-sample  distribution.  A  hypercylinder  formed  around  this  line  (with 
a  radius  typically  2  or  3  orders  of  magnitude  smaller  than  its  length  (17)) 


will  contain  a  nuaber  of  BEAST-estimated  real-sample  spectral  points.  When 
the  coordinates  of  these  points  are  transformed  into  distances  froa  the 
estimated  center  of  the  real-saaple  distribution,  a  univariate  distribution 
is  foraed.  It  is  this  univariate  distribution  that  is  used  to  construct 
confidence  liaits  by  selecting  two  quantiles  in  one  of  the  clusters  (17). 
The  reliance  of  the  BEAST  on  nonparaaetric  techniques  (techniques  that 
assuae  no  particular  underlying  distribution)  produces  an  analytical  aethod 
that  functions  without  assumptions  about  the  shape,  sise,  symmetry,  or 
orientation  of  spectral-point  clusters  in  hyperspace.  This  freedom  is 
important  because  cluster  characteristics  have  been  shown  to  be 
unpredictable  (17).  The  detection  of  abnoraal  capsules  encoapasses  both 
process -control  applications  (in  which  eapty  capsules,  rust,  aetal 
shavings,  etc.  are  of  interest)  and  the  identification  of  instances  of 
deliberate  taapering  (in  which  alaost  anything  aight  appear  in  the 
capsule).  Under  such  widely  varying  conditions  the  unpredictability  of  the 
cluster  characteristics  can  only  be  magnified,  asking  the  use  of 
nonparaaetric  techniques  even  aore  worthwhile. 


EXPERIMENTAL  SBCTION 


Equipment  Used.  Spectral  data  for  all  of  the  experiments  were  collected 
at  18  discrete  wavelengths  by  a  Technicon  InfraAlyser  400  filter 
spectrophotoaeter  connected  to  a  VAX  11/780  computer  (Digital  Equipment 
Corporation)  and  with  custom  interface,  graphics  and  database -management 
programs.  These  programs,  and  the  BEAST  algorithm  described  earlier,  were 
written  in  Speakeasy  IV  Delta  (VMS  version,  Speakeasy  Computing 
Corporation,  Chicago,  IL)  and  VAX  11  Basic  (Version  2.4).  The  programs  can 
be  obtained  by  arrangement  with  the  authors  at  the  address  given  above. 

Reproducible  positioning  of  capsules  is  important  in  minimising  the 
error  of  repeated  readings  of  the  individual  capsules.  The  initial  results 
(with  Book's  Cold  Caps  [cold  remedy  capsules])  were  obtained  by  placing  the 
capsules  into  an  elliptically-shaped  aluminum  reflector  (#1468  Progressus 
Company,  Freeport,  NY)  (see  Figure  2).  Reproducible  positioning  of  the 
capsules  within  the  reflector  was  achieved  by  removing  the  "blister*  from  a 
Cold  Caps  "blister-pack,"  trimming  it  to  about  2  mm  in  height  and  gluing  it 
into  the  center  of  the  reflector  with  the  open  side  up. 

Datril  and  Anacin-3  capsule  results  were  obtained  using  a  90°  conical 
reflector  machined  from  aluminum  (Figure  3).  A  nichrome  wire  support  was 
used  to  achieve  reproducible  upright  vertical  capsule  positioning  within 
the  cone.  Optical  surfaces  of  both  reflectors  were  polished  with  a 
commercial  polishing  paste. 

Materials  Used.  Three  brands  of  capsules  were  selected  for  this  study: 

1.  Book’s  Cold  Caps  (Book’s  Drug’s,  Inc.,  Indianapolis,  IN),  an  allergy 
and  cold  remedy  containing  a  decongestant  (phenylpropanolamine 
hydrochloride)  and  an  antihistamine  (chlorpheniramine  maleate) ,  and 
colored  red  and  white. 


2.  Extra  Strength  Datril  capsules  (500  mg  acetaminophen  [N-(4- 
Hydroxy phenyl) acetamide] ,  Bristol-Myers  Company,  New  York,  NY), 
colored  green  and  white. 

3.  Maximum  Strength  Anacin-3  capsules  (500  mg  acetaminophen,  Whitehall 
Laboratories,  Inc.,  New  York,  NY),  colored  blue  and  white. 

The  adulterants  selected  for  study  can  be  divided  into  two  categories: 
those  that  might  appear  in  capsules  as  the  result  of  process- control 
problems,  and  substances  that  are  more  likely  to  appear  as  the  result  of 
deliberate  tampering.  The  process-control  substances  tested  were: 

1.  Ferric  Oxide,  reagent  grade  (J.T.  Baker  Chemical  Company, 
Phillipsburg,  NJ) . 

2.  Aluminum  metal,  20  mesh  (Fisher  Scientific,  Fairlawn,  NJ) . 

The  substances  selected  that  are  more  likely  to  be  present  as  the  result  of 
deliberate  tampering  were: 

1.  Arsenic  trioxide,  reagent  grade  (Fisher  Scientific).  The  lowest 
lethal  dose  reported  for  humans  is  2.941  mg/kg,  or  206  mg  for  an 
average  70  kg  person  (IS). 

2.  Sodium  fluoride,  reagent  grade  (MCB,  Norwood,  OH).  The  lowest  lethal 
dose  (oral)  reported  for  humans  is  75  mg/kg,  or  5.25  g  for  a  70  kg 
person  (IS). 

3.  Crystalline  sodium  cyanide,  reagent  grade  (Aldrich  Chemical  Company, 
Milwaukee,  WI) .  The  lowest  reported  lethal  (oral)  dose  for  humans  is 
2.857  mg/kg,  or  200  mg  for  a  70  kg  person  (13). 

4.  Granular  potassium  cyanide,  reagent  grade  (Mallinckrodt,  Paris,  KY) . 
The  lowest  lethal  dose  reported  in  humans  is  2.941  mg/kg,  or  306  mg 
for  a  70  kg  person  (IS). 


All  adulterants  were  added  to  the  capsules  on  an  *as  is*  basis  without  any 
grinding,  sifting,  or  saaple  preparation.  One  or  two  grams  of  the  cyanide 
were  reaoved  froa  the  reagent  container  at  a  tiae  and  kept  covered  to 
reduce  the  absorption  of  water.  Capsules  were  filled  froa  this  covered 
reservoir  and  no  additional  steps  were  taken  to  control  water  absorption. 

Each  experiaent  performed  with  the  capsules  began  by  training  the  BEAST 
(with  10-13  unadulterated  capsules)  to  recognize  a  "good*  capsule.  The 
training  process  was  then  tested  by  using  the  BEAST  to  measure  the  distance 
(in  standard  deviations)  of  the  same  number  (10-13,  depending  upon  the 
capsule  brand)  of  "good"  capsules  from  the  center  (mean)  of  the  training- 
saaple  cluster.  The  sjurface  of  a  cluster  was  defined  as  being  3  standard 
deviations  (SDs)  away  from  its  center  (Bean) .  In  theory,  then,  all  of  the 
"good"  test  capsules  (validation  capsules)  should  appear  less  than  3  SDs 
froa  the  training-set  center,  while  all  of  the  "bad"  capsules  (contaminated 
samples)  should  appear  more  than  3  SDs  froa  the  training-set  center. 

Types  of  Experiments  Performed.  The  first  experiment  involved  placing 
arsenic  trioxide,  sodium  fluoride,  sodium  cyanide,  ferric  oxide,  and 
aluainum  shavings  in  cold  capsules  to  test  the  feasibility  of  using  NIRA  to 
analyze  the  contents  of  intact  capsules.  Two  wavelengths  were  monitored 
using  the  simple  elliptical  reflector.  On  the  basis  of  this  experiaent,  the 
second  conical  reflector  was  designed  and  tested  with  Datril  capsules. 
Further  experiments  employed  the  conical  reflector  to  analyze  intact  Datril 
and  Anacin-3  capsules  at  4  wavelengths.  The  spatial  response  and 
orientation  effects  of  sodiua  cyanide  in  Datril  and  Anacin-3  were 
deterained,  and  a  detection  limit  was  calculated  for  potassium  cyanide  in 


Anacin-3. 
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Analysis  of  the  Saaple  Reflector*.  The  empty  elliptical  reflector 
(Figure  3)  gave  log(l/R)  values  at  18  wavelengths  between  0.31  and  0.42. 
When  a  capsule  was  placed  in  the  reflector  these  values  roughly  doubled.  It 
was  thought  that  the  amount  of  specular  reflectance  passing  from  the 
reflector  to  the  detector  and  bypassing  the  sample  capsule  must  be  rather 
large.  A  new  reflector  (Figure  3)  was  constructed  to  reduce  the  specular 
reflectance.  This  60°  conical  reflector,  when  empty,  reflects  radiation 
back  toward,  its  source,  parallel  to  the  incident  beam.  When  a  saaple 
capsule  is  positioned  along  the  axis  of  rotation  of  the  conical  reflector, 
the  specular  reflectance  can  be  minimized  while  the  diffuse  reflectance  is 
maximised.  Radiation  reflected  fro*  the  surface  of  the  capsule  is  returned 
to  the  source  when  the  incident  radiation  is  perpendicular  to  the  base  of 
the  cone  (this  is  the  configuration  used  in  the  spectrophotometer) . 
Radiation  is  then  focused  along  the  length  of  the  capsule.  Any  radiation 
that  might  pass  through  the  capsules  without  being  scattered  is  also 
returned  to  the  source.  The  bulk  of  the  radiation  reaching  the  detector  is 
therefore  radiation  scattered  by  the  contents  of  the  capsule. 

If  one  assumes  that  the  base  of  the  conical  reflector  is  uniformly 
illuminated  by  collimated  radiation  (this  is  the  way  that  the 
spectrophotometer  is  designed),  the  amount  of  radiation  incident  on  any 
given  section  of  the  capsule  is  directly  proportional  to  the  curved  surface 
area  of  the  frustum  in  which  it  lies.  In  turn,  a  frustum  (a  conic  section 
taken  parallel  to  the  base  of  the  cone)  near  the  vertex  and  a  frustum  near 
the  base  of  the  same  cone  do  not  have  the  same  curved  surface  area.  (The 
curved  surface  area  of  a  frustum  is  given  by  tsCrj+rj),  where  r}  and  ra  are 
the  radii  of  the  base  and  top  of  a  right  circular  frustum,  respectively, 
and  s  is  the  length  of  the  generator  line,  i.e.,  the  length  between  the  top 
and  bottom  measured  along  the  surface  of  the  cone) . 


Suppose  that  the  length  of  the  capsule  is  divided  into  1  an  sections 
and  that  these  slices  are  numbered  fro*  1  to  20,  starting  at  the  end  of  the 
capsule  toward  the  vertex  of  the  cone.  From  the  discussion  in  the  preceding 
two  paragraphs,  the  top  slice  of  the  capsule  (i.e.,  slice  no.  20)  receives 
30  tiaes  more  light  than  the  bottom  slice  (slice  no.  1).  In  fact,  the 
aaount  of  light  (P)  received  by  a  particular  slice  numbered  R  (the  height 
of  the  section  above  the  inverted  vertex  of  the  cone)  is  given  by: 

P  =  k(2*V2)R-»V2,  (2) 

Of  course,  there  is  not  a  detector  on  each  1  on  section  of  the  capsule. 
Instead,  there  is  a  detector  inside  an  integrating  sphere,  and  the  signal 
froa  the  entire  capsule  is  integrated  to  produce  the  detector  response: 

detector  response  =  k’ (tV2)R’a-(*V2)R’ ,  (3) 

where  k  and  k’  are  proportionality  constants  that  depend  principally  on  the 
aaount  of  incident  radiation  and  the  nature  of  the  aaterial  in  the  capsule, 
and  R’  is  the  total  nuaber  of  vertical  capsule  'slices*  filled,  i.e.,  froa 
R=1  to  R’ . 

The  diaaeter  of  the  incident  beaa  in  our  instrument  is  26  aa,  asking 
direct  illumination  of  the  upper  segaents  (R=13  to  20)  by  the  incident  beaa 
the  predoainant  factor  in  producing  a  signal  froa  this  region.  The  aaount 
of  light  on  each  slice  decreases  exponentially  as  the  slice  nuaber  is 
decreased  in  this  aone.  Of  course,  the  entire  cone  is  filled  with 
scattered  light,  and  the  thickness  and  composition  of  the  capsule  wall  are 
not  uniform  over  the  capsule  length.  These  two  factors,  combined  with  the 


probable  staple  inhomogeneity ,  prevent  a  siaple  analysis  from  completely 
explaining  the  signal  observed  froa  an  individual  capsule.  However,  the 
overall  response  follows  the  trends  outlined  above. 

Coaputation  Procedure.  A  training  set  composed  of  ten  capsules  produces 
a  total  of  92378  possible  bootstrap  saaples  (calculated  froa  2n-l 

coabinations  of  n  points,  taken  n  at  a  tiae  with  replacement  froa  the 
training  set) .  Calculation  of  1000  bootstrap  replications  represents  more 
than  IX  of  the  possible  bootstrap  distribution,  a  greater  proportion  of  the 
distribution  than  is  usually  covered  by  Monte  Carlo  techniques.  A 

coaproaise  between  coverage  and  execution  tiae  Bust  be  reached  when  one 
uses  the  BEAST;  therefore,  1000  bootstrap  replications  were  used  for  all 
the  foregoing  capsule  experiments,  resulting  in  a  BEAST  analysis  tiae  of 
about  four  seconds  per  capsule.  A  BEAST  algoritha  optiaised  for  process 
control  instead  of  research  would  be  even  faster.  A  training-set  sise  of 
only  ten  capsules  is  rather  snail,  yet  is  large  enough  for  the  capsule 
experiments  because  the  variability  is  saall  froa  capsule  to  capsule  among 
the  uncontaainated  saaples.  In  general,  as  the  variability  of  the  aixtures 
in  the  training  set  increases,  the  nuaber  of  training  saaples  required  by 
the  BEAST  also  increases.  This  requireaent  is  a  common  one  aaong  NIRA 
techniques  and  is  ordinarily  not  too  burdensoae.  In  order  to  assure  that 
there  were  enough  points  (about  50  are  required)  in  the  hypercylinder  to 
set  confidence  liaits,  the  hypercylinder  radius  was  set  at  0.00060  (17). 
This  value  has  the  saae  diaensions  as  the  log(l/R)  values  collected  by  the 
spectrophotoaeter .  The  reproducibility  of  BEAST  distances  is  a  function  of 
the  nuaber  of  bootstrap  replications  employed,  the  hypercylinder  radius, 
and  the  training  set.  The  paraaeters  given  above  (and  used  to  obtain  all 
of  the  experimental  results)  have  been  found  to  give  RSD’s  around  7%  (17). 


USULTS  AND  DISCUSSION 


Tables  I  and  Table  II  suaaarise  the  results  of  the  cold-capsule  experiment. 
Ten  Hook's  Cold  Caps  were  used  to  train  the  BEAST  at  two  wavelengths  ;  it 
was  assured  that  extreme  cases  of  the  standard  capsules  were  adequately 
represented  in  the  training  set  on  the  basis  of  the  most  unusual  spectra. 
The  net  effect  of  this  procedure  is  to  make  the  training-set  cluster  in 
hyperspace  larger,  and  the  distances  measured  in  this  space  in  SDs  smaller. 
Although  this  procedure  increases  the  likelihood  that  unadulterated 
capsules  will  test  as  'good*  it  also  makes  it  more  likely  that  adulterated 
samples  will  test  similarly.  However,  this  error  would  arise  only  when 
contaminated  samples  are  spectroscopically  very  similar  to  the 
uncontaminated  ones;  such  similarity  was  not  observed  in  any  of  these 
experiments.  The  most  likely  effect  of  exaggerating  the  training  set  with 
extreme  examples  of  unadulterated  capsules  is  the  raising  of  the  detection 
limit  for  some  contaminants.  The  importance  of  this  possibility  will  be 
examined  below. 

The  distances  of  the  spectra  of  ten  unadulterated  cold  capsules,  in 
units  of  standard  deviations,  from  the  center  of  the  hyperspace  cluster  of 
training-sample  spectra  are  as  follows:  0.30,  0.69,  0.42,  0.38,  0.85,  0.70, 
0.33,  0.58,  0.23,  0.84.  These  ten  capsules  are  not  the  ten  that  were  used 
to  train  the  BEAST  and  serve  therefore  to  validate  the  results  of  the 
training  process.  The  fact  that  all  of  these  distances  are  less  than  one 
standard  deviation  indicates  that  the  BEAST  is  unlikely  to  reject  a  good 
capsule  accidentally.  Table  1  gives  the  distances  of  the  spectra  of  21 
contaminated  cold  capsules  from  the  center  of  the  training-set  spectral 
cluster.  In  every  case  the  distance  exceeds  the  3  SD  limit  set  as  the 
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dividing  line  between  "good*  and  "bad"  capsules.  The  results  indicate  that 
a  variety  of  contaainants  can  be  detected  in  intact  capsules  by  using  only 
two  NIK  wavelengths  and  a  siaple  reflector-based  capsule  aount.  Aluainus 
shavings  are  clearly  detectable  inside  the  capsule  even  though  the 
reflector  is  also  aluainua.  A  completely  eapty  capsule  can  be  easily 
differentiated  froa  a  capsule  contaainated  with  aluainua.  The  BEAST 
algorithm  works  by  detecting  the  absence  of  components  that  should  be 
present  as  well  as  by  detecting  the  presence  of  coaponents  that  should  be 
absent.  This  type  of  functioning  gives  the  BEAST  a  powerful  ability  to 
detect  all  kinds  of  taapering.  Presumably  it  is  the  absence  of  cold  reaedy, 
coabined  with  the  added  scattering  of  specular  radiation  froa  the  aluainua 
shavings  inside  the  capsule,  that  produces  the  difference  between  the 
aluainua-containing  capsules  and  the  cold  reaedy  that  is  aeasured  by  the 
BEAST. 

Spatial  Characterisation  of  the  Capsules .  When  an  adulterant  is  packed 
into  a  capsule  and  the  capsule  spectrum  is  aeasured  in  the  conical 
reflector,  the  results  can  be  fairly  well  predicted  by  our  aodel  of  the 
detector  response.  Figure  4  shows  the  actual  response  (in  SDs  froa  the 
center  of  Datril  training  set)  of  four  Datril  capsules  containing  93,  218, 
289  and  383  ag  of  sodiua  cyanide.  This  cyanide  was  packed  into  the  white 
end  of  the  capsule,  which  was  oriented  toward  the  vertex  (bottoa)  of  the 
conical  reflector.  The  slight  differences  between  the  theoretical  and 
observed  results  can  probably  be  attributed  principally  to  two  factors: 

1.  the  thickness  and  opacity  of  the  capsule  itself  are  not  constant  froa 
one  end  to  the  other  (for  instance,  there  are  two  layers  of  gelatin 
in  the  aiddle  of  the  capsule) ; 
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2.  peculiarities  are  introduced  into  the  capsule  spectra,  and  therefore 
into  the  distances  in  SDs,  by  the  particular  orientations  that  the 
capsule  contents  (both  sodiua  cyanide  and  acetaainophen  crystals) 
assuae  in  a  single  capsule  packing. 

Knowledge  of  the  characteristics  of  the  conical  reflector  peraits  soae 
spatial  profiling  of  the  capsule  for  a  given  contaainant.  A  specific  BEAST 
response  can  indicate  either  a  relatively  large  aaount  of  contaainant  in 
the  lower  sections  of  the  capsule  or  a  relatively  saall  aaount  of 
contaainant  in  the  higher  sections.  This  classic  dileaaa  -  that  of  too  aany 
unknowns  and  too  few  equations  -  is  solved  siaply  by  inverting  the  capsule 
and  running  the  BEAST  again.  The  total  aaount  of  contaainant  in  the  capsule 
is  obtained  by  developing  ordinary  NIKA  calibrations  using  hoaogeneous 
training  sets  (containing  only  one  type  of  capsule  configuration). 

Discriainating  Ability  of  the  lef lectors.  The  aost  important  feature  of 
the  reflector  cone,  however,  is  its  ability  to  discriminate  between 
capsules  with  similar  component  concentrations.  This  ability  is 
demonstrated  by  an  experiment  conducted  with  four  tainted  Datril  capsules. 
Capsules  1  and  2  contained  about  170  ag  of  NaCN,  while  capsules  3  and  4 
contained  about  460  ag  of  NaCN.  The  reaainder  of  each  capsule  was  composed 
of  the  acetaainophen  powder  normally  found  in  the  capsules,  giving  an 
average  total  capsule  aass  of  about  600  ag.  Data  were  collected  at  four 
wavelengths  with  both  the  elliptical  and  conical  reflectors,  and  the 
Euclidean  distances  aaong  the  spectra  of  the  four  capsules  was  deterained. 

The  reflector  with  the  greatest  discriainating  capability  aaxiaises  the 
average  distance  ratio 

distance  between  spectra  of  dissimilar  capsules 
distance  between  spectra  of  siailar  capsules 


This  ratio  m  calculated  for  each  reflector  using  all  possible 
combinations  of  the  four  capsules.  The  average  ratio  was  only  3.83  using 
the  elliptical  reflector,  but  diabed  to  10. S3  for  the  saae  capsules  using 
the  cone.  These  results  indicate  that  the  conical  reflector  is  sore 
sensitive  to  slight  differences  in  capsule  coaposition,  probably  because 
less  specular  reflectance  (froa  the  reflector  and  the  surface  of  the 
capsule)  reaches  the  detector  using  this  corner-reflector  configuration. 

Spectroscopically,  the  ends  of  aost  capsules  are  not  equivalent.  In 
Anacin-3  and  Datril,  the  shorter  end  of  the  gelatin  capsule  is  brightly 
colored  (blue  in  Anacin-3  and  green  in  Datril)  while  the  longer  end  is 
white  and .contains  a  light-scattering  aediua.  Figure  5  shows  the  training, 
validation,  and  NaCN-containing  saaple  sets  for  Anacin-3  in  a  two- 
wavelength  space.  The  two  distinct  training  clusters  and  the  two  distinct 
validation  clusters  are  the  result  of  including  both  of  the  possible 
orientations  of  Anacin-3  in  the  reflector  cone  (colored  end  up  and  colored 
end  down)  in  the  plot.  It  is  interesting  to  note  that  the  corresponding 
training  and  validation  clusters  are  not  necessarily  the  saae  site  or 
shape,  even  though  the  saae  capsules  were  used  in  each  and  only  the  capsule 
orientation  in  the  cone  had  changed.  In  addition,  the  aajor  axes  of  the 
NaCN-containing  clusters  are  approxiaately  perpendicular  to  those  of  the 
training  and  validation  saaples.  These  facts  deaonstrate  that  the  sise, 
shape  and  directional  orientation  of  spectral  clusters  in  space  are  not 
predictable  a  priori.  This  unpredictability,  in  turn,  violates  basic 
assuaptions  of  other  qualitative  NIRA  techniques  (£1,££),  making  their  use 
in  aany  applications  somewhat  suspect. 

Datril  capsules  containing  approxiaately  100,  200,  300  and  400  «g  NaCN 

and  in  three  configurations  (NaCN  packed  in  the  colored  end  of  the  capsule, 


in  th«  white  end,  and  nixed  into  the  acetaninophen  throughout  the  capsule) 
were  analysed  at  four  wavelengths  by  neans  of  the  conical  reflector.  It  was 
desired  to  deteraine  the  BEAST  distance  response  (in  SDs)  for  different 
aaounts  of  contaainant  in  various  locations  throughout  a  capsule.  The 
results  are  suaaarised  in  Figures  6  and  7.  The  data  in  Figure  6  were 
collected  with  the  colored  end  of  the  capsule  up  in  the  reflector  cone, 
whereas  the  data  in  Figure  7  were  collected  with  the  colored  end  of  the 
capsule  down.  In  general,  the  distances  in  SDs  neasured  froa  the  training- 
set  spectra  to  those  of  the  adulterated  capsules  are  greater  when  the 
colored  ends  of  the  capsule  are  up  in  the  reflector  cone.  The  Euclidean 
distances,  however,  are  about  the  saae  for  both  capsule  orientations.  The 
distance  in  SDs  varies  because  the  BEAST  scales  the  Euclidean  distance 
with  the  probability  of  the  point  lying  in  its  particular  direction,  and 
the  training-set  cluster  sise  (probability)  is  larger  when  capsules  are 
neasured  with  their  white  ends  pointing  up  in  the  reflector  cone.  Figure  8 
shows  one  of  six  orthogonal  views  of  the  four-dinensional  wavelength  space 
in  which  the  Datril  capsules  were  analysed.  The  Figure-8  data  were  obtained 
with  the  colored  end  of  the  capsule  up-  The  corresponding  colored-end-down 
views  are  sinilar-,  however,  the  training-set  cluster  is  slightly  larger  and 
the  distances  to  the  contaainated  saaples  are  slightly  saaller.  The  net 
effect  of  these  changes  is  to  reduce  the  sensitivity  of  NIRA  and  tbe  BEAST 
when  neasureaents  are  taken  through  the  white  end  of  the  capsule  (i.e., 
with  the  colored  end  down  in  the  cone)  ■  This  result  is  predictable  because 
NIRA  gives  inforaation  about  particle  sise  as  well  as  about  the  cheaical 
contents  of  a  capsule;  the  particle-si se  inforaation  is  obscured  soaewhat 
by  the  light-scattering  aediua  in  the  white  end  of  the  capsule. 


Figure  7  demonstrates  further  the  effect  of  taking  the  capsule-contents 
spectrum  through  two  different  layers  of  gelatin.  An  examination  of  the 
data  for  the  100  and  200  mg  samples  shows  that  packing  these  amounts  of 
NaCN  in  the  white  end  of  the  capsule  produces  greater  discrimination  (in 
SDs)  than  packing  them  in  the  colored  end.  However,  packing  300  and  400  mg 
of  NaCN  into  the  white  end  produces  a  smaller  response  than  packing  the 
same  amount  in  the  colored  end.  The  reason  that  this  reversal  is  observed 
is  simply  that  the  colored  end  allows  more  information  about  the  capsule 
contents  to  pass  through:  when  100  and  200  mg  of  NaCN  are  in  the  white  end, 
essentially  all  of  the  sample  reflectance  reaching  the  detector  passes 
through  the  white  mnd.  However,  when  300  and  400  mg  of  NaCN  are  packed  into 
the  capsule  the  capsule  is  more  than  half  full  and  a  significant  amount  of 
the  diffuse  reflectance  is  able  to  reach  the  detector  through  the  white  end 
as  well  as  the  colored  end  even  when  the  NaCN  is  packed  into  the  colored 
end  (see  Figure  0).  Unadulterated  capsules  are  quite  full  of  acetaminophen 
and  typically  some  must  be  removed  to  make  room  for  any  added  adulterant. 
The  ability  to  ascertain  the  profile  of  the  distribution  of  contaminant  in 
a  capsule  might  provide  useful  forensic  evidence  because  there  is  more  than 
one  way  to  introduce  an  adulterant  into  a  capsule.  The  distribution  profile 
has  obvious  applications  in  quality  control  as  well. 

Profiling  experiments  were  also  conducted  with  Anacin-3  capsules  using 
NaCN  as  a  contaminant.  As  with  the  Datril  capsules,  approximately  100,  200, 
300  and  400  mg  of  NaCN  were  placed  in  the  capsule  in  three  configurations 
(packed  toward  the  colored  end,  toward  the  white  end,  and  nixed  throughout 
the  capsule).  The  results  appear  in  Figures  10  and  11.  For  Anacin-3,  the 
colored  end  (blue)  of  the  capsule  might  absorb  more  in  the  near- infrared, 
relative  to  the  white  end,  than  the  colored  end  of  a  Datril  capsule 


(green).  The  ratio  of  the  scale  maxima  of  the  figures  for  Datril  (Figure 
6:Figure  7)  is  40:17  or  2.29:1,  while  the  saae  ratio  for  Anacin-3  is  20:14 
or  only  1.43:1  (Figure  10:Figure  11).  The  difference  in  these  ratios 
indicates  that  it  is  more  difficult  to  measure  the  contents  of  an  Anacin-3 
capsule  through  the  colored  end  than  it  is  to  Measure  the  contents  of  a 
Datril  capsule  through  the  colored  end.  Another  explanation  for  the  reduced 
ability  of  NIRA  to  read  through  the  colored  end  of  the  Anacin-3  capsule  is 
based  on  the  contents  of  the  capsules.  In  fact,  there  is  a  noticeable 
difference  in  the  consistency  of  unadulterated  powder  froa  a  Datril  and  an 
Anacin-3  capsule.  Anacin-3  seemed  to  consist  of  larger  flakes  than  Datril, 
and  also  had  a  greater  tendency  to  adhere  to  the  walls  of  the  capsule,  in 
spite  of  atteapts  to  empty  it.  The  amount  of  powder  remaining  in  the 
capsules  after  they  were  emptied  (but  just  before  they  were  repacked)  was 
not  measured,  and  it  is  quite  possible  that  this  aaount  was  significantly 
larger  in  the  Anacin-3  capsules  than  it  was  in  the  Datril  product.  Special 
atteapts  were  not  made  to  remove  this  clinging  powder  because  a  taaperer 
would  probably  not  make  such  atteapts  either.  'Screening*  of  the 
contaminant  by  the  acetaminophen  powder  could  therefore  be  a  significant 
factor  in  the  Anacin-3  results  observed. 

Overall,  the  histograms  show  that  the  differences  (SD  distances) 
between  the  contaminated  capsules  froa  the  training  set  are  smaller  for 
Anacin-3  than  they  are  for  Datril.  The  immediate  reason  is  that  the  Anacin- 
3  spectral  training-set  cluster  itself  is  larger  than  the  Datril  training 
set,  relative  to  the  contaminated  samples  (compare  Figure  12  for  Anacin-3 
to  the  corresponding  Datril  Figure  8).  The  fact  that  the  Anacin-3  training- 
set  cluster  is  larger  indicates  that  Anacin  3  capsules  are  normally  more 
variable  in  their  contents  than  Datrils  -  a  fact  confirmed  by  weighing  the 


capsules  in  each  training  set.  The  Bean  aass  of  the  Datril  capsules  was 

604  ag,  with  a  standard  deviation  of  5.5  ag.  The  Bean  of  the  Anacin-3 

capsules  was  670  ag  with  a  standard  deviation  of  18.2  ag.  Our  saaples  of 
Anacin-3  capsules  are  therefore  about  3.5  tiaes  more  variable  than  the 
Datril  units,  asking  the  detection  of  any  kind  of  contamination  in  Anacin-3 
a  Bore  difficult  proposition  than  the  corresponding  deteraination  in 
Datril.  Nevertheless,  the  fact  that  the  BEAST  responds  to  the  absence  of 
coaponents  that  should  be  present  as  well  as  to  the  presence  of 

contaainants  that  should  not  be  in  the  saaple  makes  the  detection  of 

adulteration  possible  under  less-than-ideal  conditions.  Table  II  gives 
BEAST  distances  in  SDs  for  NaCN-.contaainated  Datril  and  Anacin-3  capsules. 
Both  capsule  orientations  (colored  end  up  and  colored  end  down)  were 
checked  and  the  larger  of  the  two  discrimination  values  appears  in  the 
Table.  Using  the  larger  of  the  two  values  is  the  ordinary  mode  of  operation 
in  process-control  applications.  When  the  commonly  used  liait  of  3  SDs  is 
applied  to  the  training-set  cluster  it  is  apparent  that  all  of  the 
contaminated  Datril  capsules  could  be  detected  and  rejected.  All  but  two  of 
the  contaminated  Anacin-3  capsules  are  also  rejected  when  only  four 
wavelengths  are  used.  The  two  Anacin-3  capsules  that  are  not  rejected 
represent  the  lowest  NaCN  concentrations  and  in  the  aost  unfavorable 
configurations. 

The  fundaaentally  nonparaaetric  character  of  the  Quantile  BEAST  peraits 
information  vectors  other  than  near- infrared  wavelengths  to  be  used 
directly  in  the  calculations  as  though  these  vectors  were  near-infrared 
wavelengths.  For  exaaple,  the  retention  tiae  of  a  substance  in  a  liquid 
chroaatography  (1<C)  experiment  could  be  added  to  the  near- infrared 
wavelength  reflectance  data  from  n  wavelengths  to  produce  a  BEAST  analysis 


in  the  (n+1) -dimensional  space  created  by  the  addition  of  the  retention 
tiae.  Distributional  assuaptions  of  noraality  are  often  hard  enough  to 
justify  when  only  near-infrared  wavelengths  are  used,  and  the  addition  of 
dissiailar  inforaation  only  aakes  these  assumptions  wore  difficult  to 
justify.  The  perforaance  of  the  BEAST,  being  free  of  assuaptions  regarding 
data  distributions,  should  prove  to  be  even  more  superior  to  paraaetric 
aethods  in  such  applications.  The  current  proliferation  of  laboratory 
inforaation  aanageaent  systeas  aakes  a  wealth  of  inforaation  available  to 
investigators,  aost  of  which  aight  be  profitably  used  with  the  BEAST. 

In  order  for  us  to  deaonstrate  this  flexibility  of  the  BEAST  we  added 
the  total  aasses  of  the  Anacin-3  training-set  capsules  to  the  training  set. 
Capsule  aass  is  an  iaportant  paraaeter  because,  in  unadulterated  capsules, 
this  variable  is  rather  tightly  controlled.  In  addition,  the  weighing  of 
capsules  is  one  of  the  few  tests  that  can  be  perforaed  wore  rapidly  than 
NIRA.  Of  course,  capsule  weight  is  not  in  itself  a  sufficient  indicator  of 
taapering.  For  example,  unadulterated  Datril  capsules  weighed  604  ag  (SD  = 

5.5  ag)  whereas  the  NaCN-contaainated  capsules  weighed  711  ag  (SD  = 

61.5  ag) .  Accordingly,  a  substantial  portion  of  tainted  Datril  capsules 
would  pass  a  test  based  on  weight  information  alone. 

The  last  group  of  Anacin-3  distances  in  Table  II  represents  the  saae 
set  of  capsules  that  produced  the  4-D-space  distances,  except  that  the 
total  aass  of  each  capsule  was  included  to  create  a  5-D  space.  The  BEAST 
was  then  retrained  by  adding  the  total  capsule  aass  also  to  each  training 
set  saaple.  The  addition  of  the  aass  inforaation  is  enough  to  allow  the 
BEAST  to  correctly  identify  every  NaCN-containing  capsule  as  being  tainted. 

Control,  or  validation,  samples  (unadulterated  Anacin-3  capsules  that 
were  not  used  in  the  training  set)  are  also  correctly  identified  in  every 
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case  (both  with  and  without  the  weight  information)  as  being  untainted  (see 
Table  III).  Unlike  the  situation  for  Hook’s  Cold  Caps  capsules  discussed 
earlier,  the  unadulterated  Anacin-3  capsules  were  intentionally  divided 
randomly  into  training  and  validation  sets.  The  only  precaution  taken  in 
constructing  these  sets  was  to  make  sure  that  both  sets  contained 

approximately  equal  numbers  of  capsules  and  with  total  masses  below  and 

above  the  mean  capsule  mass.  The  more  random  nature  of  this  selection 
process  increases  the  distance  in  SDs  of  the  validation  capsules  from  the 
training  set.  As  in  Table  II,  the  higher  value  of  the  two  possible  capsule 
orientations  (colored  end  up  and  colored  end  down)  is  shown. 

KCN  is  perhaps  the  most  common  highly  toxic  adulterant  added  to  over- 
the-counter  drugs  (IS, IS).  The  detection  limit  for  KCN  under  optimal 
conditions  in  over-the-counter  capsules  is  thus  of  great  interest.  KCN  was 
packed  into  the  colored  end  of  eight  Anacin-3  capsules,  over  a  range  of 
concentrations  from  1  to  87%  (by  weight) .  A  strong  functional  relationship 
exists  between  the  concentration  of  KCN  in  the  capsules  and  the  distance  of 
the  capsule  (in  SDs)  from  the  training  set  determined  by  the  BEAST  (see 

Figure  13) .  This  relationship  suggests  that  the  BEAST  night  be  directly 

useful  as  a  system  or  process  control  technique  when: 

1.  the  system  is  defined  by  one  or  more  monitored  variables  (such  as  the 
wavelengths  in  this  experiment); 

2.  the  BEAST  can  be  trained  to  recognize  a  'normal*  state  as  described 
by  typical  variations  of  the  monitored  variables  (in  the  same  way 
that  the  BEAST  was  trained  by  using  a  set  of  unadulterated  capsules 
in  the  present  experiment) ; 

3.  a  given  BEAST  distance  response  in  a  particular  direction  can  be 
functionally  related  back  to  a  parameter  of  interest  in  the  system 


(as  the  quadratic  response  of  the  reflector  cone  can  be  used  to 
predict  the  location  and  amount  of  a  contaminant  in  a  capsule) . 

The  ease  with  which  the  BEAST  problem  can  be  restructured  into  a  form 
readily  solvable  by  parallel-processing  techniques  (17)  might  soon  make 
real-time  control  with  the  BEAST  an  effortlessly  attainable  goal. 

The  KCN  3-SD  detection  limit  calculated  from  the  eight  Anacin-3 
capsules  above  is  2.6  mg  (less  than  0.4X  of  the  typical  weight  of  an 
Anacin-3  capsule) .  Figure  14  depicts  the  relationship  between  the  actual 
KCN  concentration  and  the  KCN  concentration  predicted  by  NIRA  when  four 
wavelengths  are  used  for  the  eight  Anacin-3  capsules.  The  smallest  amount 
of  KCN  pladid  -in  these  capsules,  0  mg,  caused  the  capsule  in  which  it  was 
placed  to  appear  5.96  SDs  from  the  training  set  in  4-D  (4-wavelength) 
space. 

CONCLUSIONS 

The  Quantile  BEAST  method,  when  used  in  conjunction  with  NIRA  data  at 
only  four  wavelengths,  is  able  to  quickly  detect  a  wide  variety  of 
contaminants  in  capsules,  obviating  the  need  to  open  them.  The  ability  of 
this  technique  to  detect  the  absence  of  components  that  should  be  present 
as  well  as  the  presence  of  components  that  should  absent  enables  it  even  to 
signal  the  presence  of  contaminants  that  have  no  near- infrared  absorption. 
Good  results  are  achieved  on  a  simple  filter-based  instrument,  without  the 
need  for  complex  wavelength-selection  procedures.  The  detection  limit  for 
KCN  in  capsules  that  has  been  obtained  in  this  work  is  as  much  as  two 
orders  of  magnitude  below  the  lowest  reported  lethal  dose  (IS).  Substances 
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other  than  KCN  in  capsules  could  also  be  detected  at  low  concentrations. 
Selecting  analytical  wavelengths  near  the  absorption  features  of  components 
of  interest  should  improve  detection  limits  beyond  those  observed  in  these 
experiments . 

A  representative  training  set  composed  of  unadulterated  samples  is 
required  to  train  the  BEAST  algorithm  to  recognise  a  good  sample.  In  this 
research,  collecting  the  training-set  spectra  took  less  than  one  hour  and 
training  the  BEAST  algorithm  required  less  than  5  seconds.  All  of  the 
training-set  spectra  were  collected  in  a  single  day,  but  the  tampering 
experiments  took  place  over  a  period  of  two  weeks.  Nevertheless,  repeated 
runs  of  validation  samples  showed  that  the  calibration  remained  stable 
throughout  the  duration  of  the  experiment.  More  samples,  it  seems, 
collected  over  time,  would  only  enhance  the  reliability  of  the  method. 
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TABLE  I  False-Saaple  Sets  for  Hook’s  Cold  Cap  Capsules  Adulterated  with 
Various  Substances. 


Distances  froa  Unadulterated  Capsule  Training-Set  Cluster  for 
Test  Capsules  (in  standard  deviations) 


Adulterant 

Capsule  1 

Capsule  2 

Capsule  3 

Capsule  4  Capsule  5 

NaF  (100%) : 

10.26 

10.07 

10.23 

10.22  9.63 

As303  (100%) : 

10.16 

10.26 

10.33 

10.24  10.22 

A1  (100%): 

6.96 

7.24 

Al  (20%)*: 

3.58 

4.90 

Fe30,  (100%): 

8.47 

8.56 

FeaOa  (30%)*: 

6.65 

7.47 

7.12 

NaCN  (100%): 

10.57 

Eapty  capsule: 

11.09 

a.  The  reaainder  of  the 

capsule  was 

filled  with 

the  ordinary  capsule 

contents ,  i . e . ,  cold  remedy . 


TABLJ  II  Distances  of  Spectra  of  Sodiua  Cyanide-containing  Capsules  fro»  the  Training-Set 
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i»  SB.  for  «d>  c»4ul.) 


4  *»*« lengths,  4-D  .pace 


1.82 

1.38 

2.14 

1.86 
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2.38 
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1.32 
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4  wavelengths 


2.23 

2.33 

2.44 

2.74 

1.76 
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FIGUtS  CAPTIOUS 


PIG.  1.  Two  hypothetical  coapounds,  A  and  B,  whose  spectra  are  aeasured 
at  two  wavelengths.  Uncertainty  in  the  aeasureaents  at  each  wavelength  is 
represented  by  taking  1000  replicate  spectra  of  both  A  and  B,  resulting  in 
clusters  of  points  varying  about  A  and  B.  The  centers  of  the  clusters 
represent  the  best  point  estiaate  of  the  spectra  of  A  and  B.  A  line  is 
defined  by  the  centers  of  the  two  clusters,  and  the  locus  of  all  points 
within  a  user-specified  distance  of  this  line  foras  a  cylinder  in  the  space 
of  three  or  aore  diaensions. 

PIG.  2.  The  elliptical  reflector  used  in  the  initial  capsule-taapering 
experiaents.  The  reflector  fits  into  the  open  saaple  cup  supplied  with  the 
InfraAlyaer  400.  The  saaple  cup  is  then  positioned  in  the  saaple  drawer  in 
the  usual  fashion  for  analysis. 

PIG.  3.  The  80°  right-circular  conical  reflector  used  for  aost  of  the 
capsule-taapering  experiaents.  This  reflector  replaces  the  cups  provided 
for  use  with  the  spectrophotoaeter  and  fits  directly  into  the  saaple 
drawer.  The  conical  reflector  design  peraits  soae  spatial  profiling  of  the 
contents  of  the  capsule. 

PIG.  4.  The  theoretical  and  actual  response  for  a  contaainant  packed  in 
the  lower  end  of  a  Datril  capsule  in  the  conical  reflector.  The  first  bar 
((f)  represents  the  scaled  NaCN  theoretical  distance  (froa  Eq.  3)  response 
of  the  BEAST  for  an  adulterated  capsule  filled  froa  the  bottoa  to  the  top 


in  100  ag  steps .  The  second  bar  ( 1 1 1 )  represents  the  actual  response  of  the 
HKAST  for  NaCN  packed  into  the  white  (bottoa)  end  of  a  Datril  capsule. 

FIG.  5.  Spectral  clusters  of  Anaein-3  capsules  obtained  at  two 
wavelengths  with  the  conical  reflector,  and  with  both  orientations  (colored 
end  up  and  colored  end  down) .  The  cluster  on  the  left  is  foraed  by  readings 
taken  with  the  colored  end  up  in  the  cone,  while  the  cluster  on  the  right 
is  foraed  froa  the  the  saae  capsules,  read  with  the  colored  end  down  in  the 
reflector.  Training-set  capsule  (O  );  validation-set  capsule  (♦)  and 
cyanide-containing  capsule  («) . 

FIG.  6.  Datril  capsules  read  with  the  colored  end  up  in  the  reflector 
cone,  containing  approxiaately  100,  200,  300,  and  400  ag  NaCN  packed  in 
three  configurations:  in  the  colored  end  (#) ,  in  the  white  end  (|||)  aad 
aixed  throughout  the  capsule  (XXX) .  The  distance  in  standard  deviations  is 
given  in  teras  of  the  training  set  of  unadulterated  Datril  capsules  in  the 
direction  of  the  NaCN-containing  capsule. 

FIG.  7.  Datril  capsules  read  with  the  colored  end  down  in  the  reflector 
cone  (white  end  up),  containing  approxiaately  100,  200  ,  300,  and  400  ag 
NaCN  packed  in  three  configurations:  in  the  colored  end  (f) ,  in  the  white 
end  (111)  and  aixed  throughout  the  capsule  (XXX).  The  distance  in  standard 
deviations  is  given  in  terns  of  the  training  set  of  unadulterated  Datril 
capsules  in  the  direction  of  the  NaCN-containing  capsule.  The  seeaingly 
anoaalous  100  ag  'aixed*  reading  is  probably  the  result  of  particle-sise 
noise  froa,  for  instance,  a  large  NaCN  crystal  against  the  blue  end  of  the 
capsule  wall. 
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FIG.  8.  One  of  «ix  orthogonal  views  of  »  four-dimensional  space  formed  by 
taking  Oatril-capaulo  spectra  at  four  wavelengths,  for  training-set 
capsules  (♦)  and  NaCN-containing  capsules  (•) .  This  figure  shows  the 
saallest  convex  polygon  that  can  completely  surround  the  training-set 
spectral  points. 

FIG.  0.  The  capsule  orientation,  the  two  different  kinds  of  gelatin 
[scattering  (white)  and  non-scattering  (colored)]  and  the  configuration  of 
the  contents  of  the  capsule  all  affect  the  distance  response 
(discriaination  ability)  of  the  BKAST  for  a  given  contaainant 
concentration.  These  factors  can  be  used  to  advantage  to  provide  additional 
information  about  the  sample.  The  use  of  a  conical  reflector  (positioned 
with  the  base  of  the  cone  perpendicular  to  the  source  illumination  and  with 
the  vertex  down  toward  the  colored  end  of  the  capsules  in  this  figure) 
permits  spatial  profiling  of  the  capsule  for  an  identified  component. 

FIG.  10.  Anacin-3  capsules  read  with  the  colored  (blue)  end  up  in  the 
reflector  cone,  containing  approxiaately  100,  200,  300,  and  400  mg  NaCN 
packed  in  three  configurations:  in  the  colored  end  (#) ,  in  the  white  end 
(III)  and  mixed  throughout  the  capsule  (XXX) . 

FIG.  11.  Anacin-3  capsules  read  with  the  colored  end  down  in  the  reflector 
cone  (white  end  up) ,  containing  approximately  100,  200,  300,  and  400  mg 
NaCN  packed  in  three  configurations:  in  the  colored  end  (#) ,  in  the  white 
end  (III)  and  mixed  throughout  the  capsule  (XXX) ■ 


FIG.  12.  One  view  of  the  four-dimensional  space  in  which  the  Anacin-3 
capsules  were  analysed.  This  figure  corresponds  to  the  Datril  capsule 
Figure  8.  The  smallest  convex  polygon  containing  all  of  the  training-set 
capsules  is  show^  and  is  larger  than  the  corresponding  Datril  polygon.  The 
difference  in  the  sise  of  the  polygons  is  indicative  of  the  greater 
variability  of  unadulterated  Anacin-3  capsules,  a  variability  that  is  also 
reflected  in  the  weight  of  the  Anacin-3  capsules  compared  to  that  of  the 
Datril  capsules. 

FIG.  13.  The  BEAST  distance  response  (in  SDs)  as  a  function  of  the  known 
KCN  concentration  in  Anacin-3  capsules.  These  capsules  were  analysed  with 
the  colored  end  up  in  the  reflector  cone.  In  general,  the  distance  response 
is  largely  quadratic  with  concentration,  as  predicted  for  a  conical 
reflector.  The  Pearson  product -moment  correlation  coefficient  (r)  for  a 
quadratic  fit  to  the  data  is  0.085. 

FIG.  14.  Typical  analytical  working  curve  for  KCN  content  in  Anacin-3 
capsules.  The  ordinate  scale  (predicted  KCN  %)  is  found  from  a  linear 
combination  of  log(l/K)  values  at  the  four  wavelengths.  The  horisontal  line 
shows  sero  predicted  KCN  V 
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